Onset of magnetism in B2 transition metals aluminides 
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Ab initio calculation results for the electronic structure of disordered bcc Fe^^Ali-i, (0.4 < a; < 
0.75), Co2,Ali_a; and Nij^Ali-^, (i=0.4; 0.5; 0.6) alloys near the 1:1 stoichiometry, as well as of 
the ordered B2 (FeAl, CoAl, NiAl) phases with point defects are presented. The calculations were 
performed using the coherent potential approximation within the Korringa-Kohn-Rostoker method 
(KKR-CPA) for the disordered case and the tight-binding linear muffin-tin orbital (TB-LMTO) 
method for the intermetallic compounds. We studied in particular the onset of magnetism in Fe-Al 
and Co-Al systems as a function of the defect structure. We found the appearance of large local 
magnetic moments associated with the transition metal (TM) antisite defect in FeAl and CoAl 
compounds, in agreement with the experimental findings. Moreover, we found that any vacancies 
on both sublattices enhance the magnetic moments via reducing the charge transfer to a TM atom. 
Disordered Fe-Al alloys are ferromagnetically ordered for the whole range of composition studied, 
whereas Co-Al becomes magnetic only for Co concentration > 0.5. 
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I. INTRODUCTION 

The intermetallic compounds FeAl, CoAl and NiAl 
attract considerable attention due to their uncommon 
properties. For example, a high melting point makes 
them attractive as promising high-temperature aerospace 
materials. Furthermore, an unusual magnetic behavior, 
depending on temperature and concentration, has been 
observed. The /3-phase of these aluminides crystallizes in 
the B2 (CsCl) structure and persists over 45-58 at.% 
in the case of Co; /3-phases of Ni and Fe aluminides 
are stable over a broader range of composition. It is 
known that these compounds have high concentrations of 
point defects, and moreover exhibit specific defect struc- 
tures known as triple defects. EJ Predominant are appar- 
ently combinations of so-called antistructure (AS) atoms 
- usually transition metal (TM) atoms substituting Al 
(in the Al-poor region of the /3 phase) , and vacancies on 
transition- metal sites (in the Al-rich region). 

In spite of their structural similarity, these materials 
exhibit different mechanical and magnetic properties. At 
low temperatures (below 10 K), the nearly equiatomic 
alloy Coo.506Alo.494 shows a temperature dependence of 
the magnetic susceptibility which coincijies neither with 
the Curie nor with the Curie- Weiss lawE There is a com- 
mon agreement in the literature that JLhe perfectly or- 
dered compound CoAl is nonmagnetic J3^eI but already a 
small disorder leads to the onset of local moments (see, 
e.g., Ref. Is] and references therein). Also it is widely 
believed that the Co-AS atoms are responsible for the 
magnetic properties of this compound: several authors 
tried to deduce from magnetic measurements the mo- 
ment per Co-AS atom whiich is associated with an ef- 
fective moment of 5/iBBcl It is necessary to note that 
above 60 at.% Co, a permanent magnetization appears 



at low temperaturesJj For the NiAl case, no magnetic mo- 
ment has been observed for any conceaigations. Earlier 
ab initio calculations for ordered FeAE'lD reveal a mod- 
erate magnetic moment of 0.6-0.7 fi^ at the Fe site and 
a ferromagnetic ordering. Ordered CoAl and NiAl were 
found in earlier ab initio calculations to be nonmagnetic, 
which is in agreement with the measurements. 

The experimental magnetic phase diagram shows that 
FeAl near stoichiometry is in a spin-glass phase,El but in 
cold-worked disordered Fe-Al alloys the saturation mag- 
netization jsersists in both stoichiometric FeAl and Al- 
rich alloys .H 

The concentration dependence of the lattice parameter 
for the /3-phase of aluminides is also not trivial. Typi- 
cally, one should expect some slight deviations from the 
Vegard's rule, i.e. nearly linear variation with compo- 
sition. In CoAl and NiAl however, the lattice spacing 
has the maximum aJ-aoproximately equiatomic concen- 
tration of the alloy.EjO In FeAl, the lattice constant is 
nearly constant over a wide range of composition from 30 
to 51 at.% A1.EI This fact has so far no theoretical expla- 
nation. Nevertheless one may expect such behavior to be 
due to the lattice dilatation, depending on the concen- 
tration of vacancies and antisite (AS) atoms, that in its 
turn is controlled by composition and heat treatment. 

The microscopic understanding of these phenomena 
can be obtained from electronic structure calculations. 
Actually the electronic structure of magnetic TM alu- 
minides with stoichiometric composition has been al- 
ready studied for many times by various methods in 
different approximations. In the earliest -calculations 
within the KorringaJiohn-Rostoker (KKR)E3 and mod- 
ified KKR methodsE^I the problem of filling up of the 
TM d-bands by Al p-electrons was discussed in detail 
and the charge transfer was shown from an Al to a TM 



site. The possibility of magnetic ordering was studied 
for the Ni-Al system by the augmented spherical wave 
(ASW) methodtil. Hereby only for NisAl a ferromag- 
netic ordering was found, with a value of the Ni mag- 
netic moment by an order of magnitude less than in pure 
Ni. Optical properties of TM aluminides were studied by 
the linear mufRn-tin orbital (LMTO) methodta and bx 
the linearized-augmented-plane-wave (LAPW) methoclE£l 
for NiAl and CoAl. The experimental absorption max- 
ima were fouad to be correlated to the band structure of 
these alloys.llZl Further details of the electronic structure 
of these compounds were studied by comparing the data 
from electron energy-loss spectroscopy with theoretical 
investigations performed by the LMTO method.li3 The 
trends in the chemical bonding and the phase stability 
of transition metal aluminides with equiatomic compo- 
sition have been studied by the full-potential linearized 
augmented plane-wave (FLAPW) methodE A review of 
electronic structure calculation results, along with band 
structures, densities of states and Fermi surfaces of many 
TM aluminides can be found in Ref.Ej. A recent study 
using the full-potential linearized augmented Slater-type 
orbital methodc3 reports the formation energies and equi- 
librium volume of many 3d aluminides. 

The equation of state and all zerctpressure elastic 
moduli for CoAl have been calculatedE^ also using the 
FLAPW method. In addition, cohesive, electronic 
and magnetic properties of the transition metal alu- 
minides have been calculated using the tightrbinding lin- 
ear mufRn-tin orbital (TB-LMTO) method.El Here it was 
found that only FeAl retains a magnetic moment, which 
has been mentioned above. These findings coincide with 
the results of earlier LMTO jmlculations for NiAl and 
FeAl intermetallic comppundsEj and with ab initio pseu- 
dopotential calculationsE3 for CoAl. 

There are relatively few calculations aimed at the 
study of defects' influence on the electronic and mag- 
netic structure of TM aluminides. Ab initio electronic 
structure calculations for point defects in CoAl have 
been performed by Stefanou et alB within the KKR- 
Green's function method. Herein the perturbation of 
the potential was included at the sites which are near^ 
to the substitutional impurity. Earlier, Koenig et al.cJ 
performed Green's function calculations for vacancies in 
B2 aluminides by the LMTO-Green's function method, 
but allowing only the potential at the impurity site to 
be perturbed. The LMTO method has been applied 
to study the electronic structure of AS defects in FeAl 
where the point defect was modelled by suitably cho- 
sen supercellsrj Finally, the LMTO-CPA technique has 
been used for the calculation of electronic and_thcrmo- 
dynamic properties of disordered NiAl alloys. E3 In Ref. 
P7| , the LMTO-CPA calculation results are presented for 
discussing the order-disorder transition in FeAl alloys. 

The supercell approach has been used in ordjer.to study 
the antiphase boundary in NiAl and-.FeAlL3n3 as well 
as point defects in these aluminides.E3 Very interesting 
results for the magnetism of Co-doped NiAl by FLAPW 



supercell calculation have been obtained by Singh,c2l and 
we compare our findings with these results. 

Although many studies have been performed up to 
now, the aspects of magnetism seen as a function of con- 
centration remains not completely clear in these com- 
pounds. The purpose of this paper is therefore to inves- 
tigate the onset of magnetism in the /3-phase of ferro- 
magnetic TM (Fe, Co, Ni) aluminides. First-principles 
total energy calculations have been performed in order 
to determine the magnetic phase stability of these al- 
loys. The calculations are based on the density-functional 
(DF) theory. We have applied the TB-LMTO method 
in order to solve the DF equations in the case of near- 
stoichiometric ordered compounds with points defects. 
The KKR-CPA has been used for disordered substitu- 
tional alloys. We also present here the results of our su- 
percell simulations for AS atoms and vacancies for both 
sublattices of the ordered compounds and compare them 
with the results for completely disordered alloys in order 
to clarify the physical reasons for the magnetic instability 
in these systems. 

This paper is organized as follows: Sec. II is devoted 
to the description of our calculational methods, in Sec. 
Ill we discuss the results of our TB-LMTO calculations 
for the point defects. In Sec. IV we present the results 
for disordered alloys, concentrating on Co-Al as the most 
interesting one. A summary is given in Sec. V. 



II. METHODS OF CALCULATION 

The conventional linear muffin-tin orbital method in 
the atomic spheres approximation (LMTO-ASA) as well 
as its tight-binding formulation_CTB-LMTO) have been 
well described in the liter ature.E2lc3 Here we present only 
the relevant computational details. 

We have used equal and space-filling atomic sphere 
radii on Al and TM sites in all calculations, scaling 
them with the lattice constant when studying the volume 
trends. According to the experience of previous studies 
(see, e.g. Ref. ^), this choice gives the best description 
of electronic densities for these systems in the ASA. Basis 
functions up to ^ = 2 were included for both constituent 
atoms explicitly, aad the I — 3 functions via the down- 
folding procedure. a Combined corrections have been in- 
cluded throughout, and the effect of non-local corrections 
(after Langreth and Mehl, Ref. ^5|) to the local density 
approximation (LDA), i.e. compared with the-exchange- 
correlation potential by von Barth and Hedin,E3 has been 
specially addressed and discussed. Since the advent of 
schemes incorporating the non-local functional forms of 
the exchange-correlation energy in terms of electron den- 
sity, they have been applied to many solid-state prob- 
lems, including magnetic ones. Whereas the non-local 
corrections were shown to improve the LDA-based de- 
scription in many cases (providing, e.g., the correct to- 
tal energy hierarchies of magnetic phases), one should 



keep in mind that they tend to underestimate the bond- 
ing and hence to lead to an overestimated equihbrium 
volume (otherwise typically underestimated within the 
LDA). Moreover, the non-local corrections favor larger 
magnetic polarization. In some systems, where the in- 
terplay of magnetic properties and volume is crucial, 
one must be careful in applying the non-local-corrected 
schemes. In the present study, we in many cases refer 
to the results obtained both within the LDA and incor- 
porating the non-local corrections. Taken together, they 
provide reasonable error bars for the properties obtain- 
able from a first-principle calculation spheme. 

The tetrahedron method integrationEZI have been per- 
formed on the 16x16x16 mesh over the Brillouin zone 
(BZ), so that the total energy values cited below were 
converged against the fc-mesh enhancement within rele- 
vant digits, important for the study of the volume trends. 

The electronic structure of disordered Co-Al al- 
loys havCj— been calculated using the fast KKR-CPA 
techniqueEa The KKR-CPA method is based on the mul- 
tiple scattering theory. In this case the electronic Green's 
function has to be calculated self-consistently with the 
scattering path operator for an atom embedded into the 
ordered lattice of effective scatterers. In addition, the 
potential for each constituent has to be brought to self- 
consistency in the CPA sense. Therefore in KKR-CPA we 
have two self-consistent procedures, that leads to a num- 
ber of numerical difficulties. An effective solution of this 
problem has been proposed in Ref. Bq. Specifically, for 
the BZ integration the tetrahedron scheme with complex 
energies is applied, and in the CPA loop the extrapolation 
approach with only few iterations is introduced in order 
to solve the CPA equation on the rectangular contour 
in the complex energy plane. As a rule, we used in the 
present work the same parameters for the energy integra- 
tion contour as described in Ref. ^. The necessarily min- 
imal basis set included s, p, d and / orbitals for both TM 
and Al atoms, at which all space- filling atomic spheres 
had the same size. In this sense, the calculation setup 
was close to that used in TB-LMTO calculations for the 
ordered structures. Another similarity is that no frozen 
core approximation has been used, i.e. both schemes are 
all-electron ones. The difference is, TB-LMTO calcula- 
tions were scalar-relativistic ones whereas KKR-CPA cal- 
culations were nonrelativistic, but this difference seems 
to have only negligible effect on structural and magnetic 
properties of TM aluminides. Also, the effects of non- 
local density functionals have not been studied in our 
KKR-CPA calculations. The exchange-correlation has 
been treated there in the LDA according to the Hedin- 
Lundqvist prescriptionEa for non-magnetic SKStems and 
with spin scaling after von Barth and HedinEEl for mag- 
netic alloys. 

In contrast to the LMTO method where all energy 
eigenvalues are explicitly provided by matrix diagonal- 
ization and the application of the tetrahedron method 
is straightforward, the BZ integration is more involved 
in the KKR method. In order to trace the a priori un- 
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FIG. 1. Local DOS of components in the non-magnetic case 
for ordered FeAl, CoAl and NiAI as calculated by TB-LMTO. 
The dashed line indicates the Fermi level. 



known number of poles in the energy-dependent determi- 
nant, we have divided the irreducible part of the BZ into 
1024 equivalent tetrahedra by an uniformly distributed 
mesh of /c-points. While examining the tetrahedra one 
by one, we have divided those where the phase of the in- 
tegrand changes by more than 7r/3 over the volume into 
eight smaller tetrahedra, etc. This subdivision allows us 
to reduce the number of fc-points necessary to obtain the 
relevant precision by a factor of 2-2.5. 

In consequence, 10-20 iterations using the modified 
Broyden schemec2l for accelerating the convergence have 
been sufficient in the potential loop, at which the CPA 
self-consistency being achieved inside each iteration for 
the potential. The charge self-consistency has been used 
as the convergency criterion, and the iterations have 
been performed until the electronic density differences 
were stable to within 10^^. Apart from conventional to- 
tal energy and electronic structure calculations within 
the CPA,|-we have used the fixed-spin-moment (FSM) 
procedurecll in order to analyze the volume/magnetic mo- 
ment total energy dependence in Co-Al alloys. 



III. ELECTRONIC STRUCTURE OF POINT 
DEFECTS IN CoAl AND FeAl COMPOUNDS 

The above cited previous calculations of ordered FeAl, 
CoAl and NiAI compounds have shown that their elec- 
tronic structures are rather similar and related roughly 



by a rigid-band shift, which is depending on the number 
of conduction electrons and hence on the placement of the 
Fermi energy. This is illustrated by the electronic densi- 
ties of states (DOS) for these materials, which are shown 
in Fig. n^. The fcrromagnetism of FeAl is perfectly consis- 
tent with the Stoner model, as the Fermi level crosses the 
highest peak of a pronouncedly two-peak DOS distribu- 
tion. For CoAl and NiAl, the Fermi level is outside of the 
region which is favorable for the onset of ferromagnetism. 

Singho emphasized the existence of a peak in the 
"rigid" DOS between the positions of the Fermi level in 
CoAl and NiAl. One can expect that tuning the num- 
ber of conduction electrons, e.g. doping NiAl with Co, 
would favor the appearance of the magnetic order. In 
Ref. ^, for example, the effect of substitutional Co in 
NiAl has been investigated in a sequence of supercell 
calculations with the FLAPW method. An instability 
against ferromagnetism was found for Coo.25Nio.75Al and 
Coo.i25Nio.875Al compositions. Apart from this, strong 
deviations from the rigid-band behavior were found for 
doped systems, implicating that the onset of magnetism 
in these systems may also be strongly dependent on the 
short-range order. 

In the present study, we are primarily interested in the 
influence of local defects (vacancies and antisites). 



FeAl 
Fe antisite 



CoAl 
Co antisite 




-0.8 -0.4 0.0 0.4 -0.8 -0.4 0.0 0.4 
Energy (Ry) 

FIG. 2. Local DOS at the AS Fe defect and its several 
neighbors in FeAl as calculated by TB-LMTO. The dashed 
line indicates the Fermi level. 
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FIG. 3. Local DOS at the AS Co defect and its several 
neighbors in CoAI as calculated by TB-LMTO. The dashed 
line indicates the Fermi level. 

which we have introduced on both sublattices in FeAl 
as well as in CoAl. In order to suppress the defect- 
defect interactions, quite large supercells have been used. 
For examining the energetics (i.e., optimizing the lattice 
spacing for different defect configurations), we have used 
supercells with translation vectors (0, 2a, 2a), (2a, 0, 2a), 
(2a,2a,0) of the underlying cubic lattice, i.e. including 
32 atoms. For reasons of controlling and for the thor- 
ough study of the spatial charge and magnetic moment 
distribution, we have performed several calculations with 
twice larger supercells, spanned up by the translation 
vectors (—2a, 2a, 2a), (2a, —2a, 2a), (2a, 2a,— 2a). In the 
last case, the next translated defect appears only as the 
9th nearest neighbor within the corresponding sublat- 
tice. Since the calculations are fully self-consistent, they 
do probably provide a more accurate description of the 
screening of defects as was previously ppbtained with the 
use of KKR-Green's function method,El where the poten- 
tials only at the defect site and its nearest neighbors were 
allowed to deviate from the perfect-crystal values. 

It should be noted that in the present treatment we did 
not account for a possible lattice relaxation of neighbor- 
ing atoms around a vacancy of any type, keeping in mind 
that such relaxation cannot be reliably addressed in the 
ASA. To our knowledge, no such relaxation studies (with 



TABLE I. Uncompensated number of electrons in the atomic sphere Q and local magnetic moments m in perfect ordered 
FeAl, CoAI and over shells of atoms around point defects in these alloys. The results are from 32-atom supercell calculations 
with nonlocal exchange-correlation potential and a=5.40 a.u. 
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the use of any full-potential scheme) have yet been done 
for the systems in question. According to Ref. ^, the 
relaxation around Pd vacancies in PdAl is up to 2.5%, 
and expectedly even less in FeAl, i.e. too small to have 
a noticeable effect on magnetic moments. 

As most interesting examples related to the develop- 
ment of magnetic properties, the local DOS at different 
spheres around Fe and Co AS impurities in FeAl and 
CoAl are shown in Fig. g and y, correspondingly. The 
numerical data for charge transfer values and magnetic 
moments inside atomic spheres of all defect systems stud- 
ied are presented in Table 0. These data correspond to 
the lattice constant a=5.40 a.u. that is about the equi- 
librium value for both ordered FeAl and CoAl. In the 
discussion that concerns the bulk properties and defect 
formation energies (see below) we refer to optimized equi- 
librium lattice spacings for defect systems, but this is not 
so important for discussing partial charges. The asymp- 
totic value of the charge transfer within the perfect bulk 
is roughly recovered in the crystal away from defects. 
Some fluctuations in both charge transfer and magnetic 
moments however prevail and are discussed below. 

It is well known that the charge transfer is not uniquely 
determined, since it depends on the choice of atomic 
spheres or cells between which the charge flow occurs. 
Moreover, because of different localization degree and 
spatial distribution of charge in different solids it doesn't 
make sense to discuss the charge transfer accompanying, 
e.g., the formation of a compound from elemental con- 
stituents. This has been discussed at length by Schulz 
and Davenport for "id aluminides in Ref. ^. With this in 
mind, we compare in Table | the uncompensated charges 
Q in equal-sized atomic spheres at Al and TM sites in 



perfectly ordered aluminides. Charge- induced perturba- 
tions (Friedel oscillations) around the defects also can be 
found in Table |. We emphasize that not the absolute 
values of Q, but their changes over shells of neighbors to 
a defect are meaningful and of primary interest here. The 
interesting feature seen in Table | is that both Fe- AS and 
Al-AS defects are less charged as compared to respective 
bulk atoms of the ordered alloy. This is not surprising 
because the electro-negativity of the nearest neighbors is 
the same in both cases. 

The general trends for the Fe-Al system could be al- 
ready expected based on the rigid-band model. Since the 
electron transfer in the ordered alloy occurs from Al to 
TM, the ways to lower the concentration of electrons at 
the TM site are substitution of Al with a TM atom, or 
vacancies of any kind. The Al-AS substitution, on the 
other hand, is unfavorable for magnetism, as can be seen 
from Table 11. The Fermi energy in ordered FeAl crosses 
the upper slope of a high peak. According to the rigid 
band picture, the decrease in the Al— i-Fe charge trans- 
fer leads to an increase of the local DOS value at the 
Fermi energy and thus favor the increase in the magnetic 
moment in the framework of Stoner theory. 

Consistently with this it is seen from Table | that the 
magnetic moment of Fe as an AS defect (see also Fig. 0) 
or as a neighbor to the Al vacancy is enhanced simultane- 
ously with the decrease of the charge transfer onto these 
sites. Because of such a locally changed charge transfer 
the local DOS at the defect site is affected in a rather 
complicated way, different from what one would expect 
from rigid-band considerations. The local DOS at the AS 
Fe, which is the central atom of the Feg nearest-neighbor 
cluster, acquires certain similarity with the DOS of pure 



bcc Fe with its half-filled minority-spin subband. Fur- 
thermore the magnitude of the local magnetic moment is 
quite close to that of bulk Fe. This is consistent with the 
fact that practically no charge transfer occurs between 
the central Fe atom pad its nearest neighbors. From the 
experimental studiescj the local moment has been esti- 
mated to fall within the range of 4 to 5.4 /is, depending 
on the vacancy concentration. According to our calcula- 
tion, the cluster of nine Fe atoms including the AS de- 
velops 8.3 Hb- The local DOS of Fe atoms which are 
more distant from the defect [Fe(iii) and Fe(|||) in 
Fig. H become rapidly resembling that of Fe in ordered 
magnetic FeAl. It is noteworthy that the spin-up sub- 
band is essentially filled, so that any subsequent removal 
of electrons may take place primarily from the minority- 
spin subband and thus increase the magnetic moment. 
Therefore one should expect a clear correlation between 
fluctuations of charge transfer and those of local mag- 
netic moment over shells of distant Fe neighbors. Apart 
from the Fe-AS atom, the pronounced enhancement of 
the magnetic moment, compared to the average value 
in the ordered alloy, occurs for Fe atoms which are first 
neighbors to the vacancy at the Al site. 

From all four defects considered in the Co-Al system, 
only the Co-AS develops a magnetic moment which is 
in agreement with previous observations based on the. 
KKR-Green's function calculations of Stefanou et alM 
Quantitatively, our results differ only slightly from those 
of Ref. M our magnetic moment at the AS Co (1.53 /is) is 
almost the same as in pure elemental Co (1.72 fi^ exper- 
imentally; 1.55 Ub as calculated, e.g., in Ref. pi) against 
1.22 /Xb in Ref. g. The total moment at the cluster of nine 
Co atoms is 2.67 /Xb (2.06 /Xb in Ref. |5|). This difference is 
probably related to the fact that the magnetic moments 
are attributed to space-filling atomic spheres, which we 
used in our case, compared to non-overlapping muffin-tin 
spheres in the KKR-Green's function approach. 

An unusual feature in Table | is that in two cases, 
namely for the Fe vacancies and for the Al AS in FeAl, 
the Fe atoms which are most distant from the defect ex- 
hibit quite large magnetic moments. This is most prob- 
ably an artefact related to a yet too small supercell size. 
This artefact, which is quite stable against the changes 
in the calculation setup (more dense fc-mesh in the BZ 
integration etc.), was not reported in e|aElier LMTO su- 
percell calculations by Gu and Fritsche.E3 They used the 
same geometry of the 32-atoms supercell (for AS defects) 
as we did and of course obtained much similar results for 
the magnetic moments in the vicinity of defects. Our 
explanation for the artefact is the following: Magnetic 
moments within different shells of Fe neighbors surround- 
ing the defect exhibit long-range oscillations around the 
bulk value of 0.76 /Xb, as is seen from Table |. For the 
Fe atoms most distant from the vacancy at Fe and from 
the Al AS, which are both situated at 2a from the defect, 
the oscillations just tend to be positive. In our partic- 
ular supercell geometry there is only one of such an Fe 
atom per supercell, which at the same time is the (fourth 
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FIG. 4. Charge transfer (left scales) and magnetic mo- 
ments (right scales) over shells of Fe neighbors to AS Al and 
Fe vacancy in FeAl. Dots - calculation results for 32-at. su- 
percells; open circles - for 64-at. supercells. Dashed lines 
indicate the asymptotic values of charge transfer and mag- 
netic moment in ordered FeAl. 



nearest) neighbor to six defect sites. The constructive 
interference of these oscillations coming from six defects 
is pinned at a single atom, resulting in an anomalously 
high magnetic moment. In a real defect crystal with no 
ordering in the mutual orientation of defects there is no 
physical reason for such a high moment to be pinned at 
any site far from the defect. 

In order to prove that the explanation given above is 
correct, we have performed the calculations for super- 
cells of increased size (including 64 atoms) for two defect 
types in question. The magnetic and charge distributions 
over spheres are shown in Fig. in comparison with the 
data for the 32-atom supercells, which are listed in Ta- 
ble |. One can see that, while the anomalous magnetic 
moment at the fourth neighboring Fe site is considerably 
reduced in a larger supercell, the trends in the vicinity of 
defects are not as much affected. Therefore the use of 32 
atom-supercells for our subsequent analysis is justified. 
An additional observation seen from Fig. ^ is that charge 
fluctuations and magnetic fluctuations over spheres be- 
have almost identically. This is due to the manifestation 
of the above mentioned feature that local charge fluctu- 
ations affect essentially the minority-spin subband and 
hence the net magnetic moment. 

For CoAl, the Fermi energy is placed near the mini- 
mum between two d-peaks and therefore leads to a loss of 
magnetic moment in spite of decreasing of charge transfer 
due to the charge-induced perturbations around the Co- 
vacancy, Al-vacancy and Al-AS defects. Only the Co-AS 
defect provides a sufficient decrease in the charge transfer 
in order to fulfill the Stoner condition and to create the 
magnetic moment at this site. 

It is necessary to point out that planar defects, like 
the surface, may be expected to affect the magnetism in 
a similar way as point defects do. In particular, FLAPW 
7-layer slab calculations, performed in ordeiLto simulate 
the surfaces of FeAl and CoAl compounds ,l3 resulted in 



TABLE II. Cohesive properties (eV) at equilibrium lattice constants a (a.u.) of FeAl, CoAI and NiAI 







FeAl 




CoAl 




NiAI 






Energy 




a 


Energy a 


Energy 
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Experiment (Ref. [45-] 
and as cited m Rcfrg) 

TB-LMTO (Ref. |) 


-0.26 
-0.50 




5.409 
5.364 


Formation energy 
-0.56 5.408 
-0.75 5.317 


-0.64 
-0.77 




5.456 
5.377 


Present work: Ordered B2 


-0.40 




5.397 


-0.66 5.354 


-0.71 




5.422 


Present work: Disordered 


-0.23 




5.366 


-0.38 5.367 


-0.53 




5.410 


Ordered B2 


-0.02 






Magnetic energy 








Disordered 


-0.05 






0.00 
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a magnetic ordering for FeAl(lOO) and CoAl(lOO) sur- 
faces, while NiAl(lOO) was magnetically dead. The Fe 
surface atom with a magnetic moment of 2.57 /is has up 
to 0.25 electrons less than the next iron layer atom with 
a moment of 0.67 /is - this is the same trend, but more 
pronounced, as in our calculation for a Fe atom neighbor- 
ing the "cavity" of the Al vacancy. In the case of CoAl, 
the surface atom loses about 0.3 electrons as compared 
to the atom in the next Co layer, which results in the 
enhancement of the surface Co moment to 1.12 fi^. 

The changes in the total energies due to the variations 
of the lattice constants exhibit the expected parabolic be- 
havior. This is shown in Fig. ^, where the results for two 
different exchange-correlation potentials are presented. 
Table || includes the results obtained by a parabolic fit- 
ting of the total energy data. As pointed out in Sec. 0, 
the data related to ordered compounds, as well as the 
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FIG. 5. Total energy (dots, left scale) per formula unit vs. 
lattice constant in ordered FeAl and CoAl as calculated by 
TB-LMTO in the LD A and with non-local corrections (LMH) . 



defect formation energies, have been obtained with non- 
local corrections to the exchange-correlation, whereas 
the data for disordered alloys (discussed in more de- 
tail below) result from the calculations within the LDA. 
The calculated equilibrium lattice constants are consis- 
tently underestimated by about 1%, with somehow better 
agreement to the experimental data in the case where the 
non-local corrections have been included. The compaij^ 
ison with the results of earlier TB-LMTO calculations^ 
for ordered TM aluminides, obtained within the same 
formalism but in the LDA, reveals that the non-local 
corrections to the LDA produce indeed a systematic im- 
provement of all cohesive properties for the alloys under 
consideration. 

The cohesive properties have been calculated at equi- 
librium volumes and are compared with the experimental 
data available. The experimental data in Table || are 
from thermochemical formation energy measurements; 
the total energy of the constituent solids are taken at 
the equilibrium volume in their respective ground-state 
structures, i.e. fee Al and Ni, hep Co and bcc Fe (ferro- 
magnetic Fe, Co and Ni), calculated with the same setup 
as used for alloys. The total energies in the TB-LMTO 
calculation have been also minimized for the defect super- 
cell systems, thus accounting for a uniform lattice dilata- 
tion around defects. The corresponding lattice constants 



are shown in Table HI 



It is noteworthy to mention that the calculated heats 
of formation come out systematically overestimated for 
all ordered compounds, but to a less extent iruthe present 
work than in earlier TB-LMTO calculationsO with local 
exchange-correlation. On the other side, they are sys- 
tematically underestimated by about the same margin in 
the KKR-CPA calculations (see next Section). Appar- 
ently, the properties of real systems fall between the two 
extremities of complete substitutional disorder and com- 



TABLE III. Equilibrium lattice constants a (a.u.) from 
TB-LMTO supercell calculations 

TM vacancy Al vacancy TM antisite Al antisite 



FeAl 
CoAl 



5.404 
5.355 



5.401 
5.337 



5.387 
5.341 



5.424 
5.381 



plete ordering with point defects, which are addressed in 
the present study. 

In our calculations ordered CoAl and NiAl fail to de- 
velop any magnetic ordered state in agreement with all 
previous calculations, while FeAl has both ferromagnetic. 
(FM) and nonmagnetic solutions. The earlier LMT0E3 
and TB-LMTOQ studies found the nonmagnetic solution 
as the ground state within the LDA; our calculation with 
the LDA puts the magnetic solution minutely lower in en- 
ergy, actually one cannot reliably claim this (see Fig. ||). 
With non-local corrections included to the exchange- 
correlation, the ferromagnetic ground state clearly gains 
in energy with respect to the nonmagnetic state. The 
values of the magnetic moment (per cell) for different 
values of lattice constant are also shown in Fig. |^. As is 
consistent with other studies on magnetic systems, the 
non-local exchange-correlation results in somehow larger 
values of magnetic moments. It shouldJpe noted also that 
ASW fixed-spin-moment calculationg^ provided, in ad- 
dition, the antiferromagnetic (AFM) state with an energy 
difference between FM and AFM states of only a tenths 
of a mRy. Since three possible ground states have very 
small energy differences, one can conclude that the for- 
mation of a spin-glass state seems plausible in the Fe-Al 
system. The existence of the spin-glass state was indeed 
detected in Ref. ^. 



IV. ELECTRONIC STRUCTURE OF 
DISORDERED ALLOYS 

While the supercell approach aims at studying point 
defects in otherwise perfect alloys, the limiting case of 
the ultimate substitutional disorder can be treated by 
the KKR-CPA approach. Here an atom of each species 
in the lattice with the probability equal to its concen- 
tration can be considered as an impurity in the effec- 
tive medium, and with complementary probability - as 
the AS defect. We have performed such calculations for 
Fea;Ali_a;, Coa:Ali_a; and NiajAli-^, alloys in the range of 
compositions x = 0.4 to 0.6. From the series of KKR- 
CPA calculations for different alloy compositions and dif- 
ferent volumes per atom we have determined the equilib- 
rium lattice parameters. Their values for x — 0.5 are 
also listed in Table ||. In Fig. ^ we show the results of 
these calculations for Co-Al and Fe-Al systems for both 
ferromagnetic and nonmagnetic solutions. The magnetic 
solution is missing in Co^^Ali-a; only for x — 0.4. For 
Nij^Ali-^, no magnetic solution has been found in the 
whole range of concentrations, which have been studied 
in this contribution. 

As it is shown in Fig. ||the C00.5AI0.5 alloy exhibits the 
ferromagnetic solution, in contrast to the ordered CoAl 
compound (see Fig. ^. The total energy difference be- 
tween ferromagnetic and nonmagnetic ground states is 
however extremely small, as is listed in Table ||. The 
magnetic moment at the Co site is 0.47 /ig near the tran- 
sition point to the nonmagnetic state. The ferromagnetic 
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FIG. 6. Total energy (per atom) vs. lattice constant in 
disordered Fe-Al and Co-AI for several concentrations as cal- 
culated by KKR-CPA, for paramagnetic (PM) and ferromag- 
netic (FM) solutions. 



state disappears under pressure but becomes clearly the 
favorable one for increased volumes in this system. 

In order to investigate the type of magnetic instability, 
we have performed the FSM analysis for the C00.5AI0.5 
alloy, with the total energy calculated as a function of 
fixed magnetic moment and lattice constant. The re- 
sults of these calculations shown in Fig. [^ indicate the 
presence of a saddle point between magnetic and non- 
magnetic solutions. It means that the phase transition 
between ferromagnetic and nonmagnetic states occurs as 
a first-order transition. At the moment, this is only a 
theoretical prediction, since the fully disordered CoAl al- 
loy was not obtained experimentally and, moreover, the 
magnetic transition is expected to occur for a somehow 
expanded lattice. However, this transition can be seen 
as a kind of pressure-induced transition in the case of Co- 
rich disordered alloys, as follows from Fig. |^, and it has 
been jeally found at low temperatures in the Cog eAlo 4 
alloyi 

For Fe-Al disordered alloys, the ferromagnetic solu- 
tions exist over the whole range of composition 0.4 < 
X < 0.75, which has been covered in the present work. 
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TABLE IV. Disordered alloys formation energies (eV) 
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FIG. 7. Total energy vs. lattice constant and magnetic mo- 
ment calculated by FSM KKR-CPA for disordered C00.5AI0.5. 

The results are shown in Fig. |6[ The energy difference be- 
tween ferromagnetic and nonmagnetic solutions increases 
for Fe-rich disordered alloys, and for the Feo.sAlo.s alloy 
we obtain the ferromagnetic ground state already within 
the LDA approach (see, also. Table [V). Only for the case 
X < 0.4 the nonmagnetic solution can be more stable. It 
is interesting to point out that for both Co-Al and Fe-Al 
disordered alloys the minimum of the formation energy 
coincides with the stoichiometric composition. On the 
contrary, for the Ni-Al alloy the minimum is shifted to a 
higher concentration of the transition metal. 

The Fe-Al alloys, annealed from high temperatures in 
order to obtain the perfect disordered state, show the 
magnetization in a broad range of compositions.^ We 
compare these experimental values with the calculated 
KKR-CPA magnetization per unit cell in Fig. 0, where 
also the local magnetic moments for Fe and Al atoms are 
given. Actually, the calculated magnetization decreases 
slower with the Al concentration than it does in the ex- 
perimental findings. This can be explained by the change 
of disordering degree with concentration in the actual 
samples. The comparison between the experimental and 
the calculated dependencies of the lattice constant sup- 
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ports this possibility (see Fig. ||) . The experimental de- 
pendence of the lattice constant on composition can be 
considered as an averaging over local variations of spac- 
ing around randomly distributed defects, as affected, e.g., 
by the heat treatment. The KKR-CPA calculation, on 
the contrary, gives the equilibrium volume as a smooth 
monotonic function of the iron concentration. As usual, 
the magnetization is a function of the equilibrium vol- 
ume, and an extended lattice provides a larger value of 
the magnetization. In order to get more insight into the 
volume trends on disordering, the data on equilibrium 
lattice constants for defect systems, as estimated from 
our supercell calculations (see Table [II), may be of some 
use. We did not allow for local (breathing) relaxations of 
neighbors around the defect because such energy evalu- 
ations might be not sufhciently reliable in a calculation 
done with the ASA. The trends in the uniform lattice ex- 
pansion, simulating in this case the effect of a lattice di- 
latation around defects, are usually quite robust and rea- 
sonably obtainable in the ASA. The equilibrium lattice 
constant for a system with an Al-AS defect is 5.424 a.u., 
whereas that with an Fe-AS defect is 5.387 a.u. The only 
possibility to understand the experimental trend seen in 
Fig. H, i.e. that the lattice constant is relatively indepen- 
dent on the Fe concentration, - is to assume the existence 
of Al antisite defects even at elevated Fe concentration. 

It is very useful to compare the DOS for the disor- 
dered alloys series Ni-Al, Co-Al and Fe-Al as we do 
in Fig. y for ordered compounds. This comparison is 
shown in Fig. pi As becomes evident from our results, 
the rigid-band model applies better for any fixed TM-Al 
composition with different TM components than for the 
same alloy with different concentrations. In particular, 
the low-energy structure at about —0.3 eV is present for 
Ni-Al, Co-Al and Fe-Al alloys with 60% of TM but it 
disappears for the alloys with smaller TM concentrations. 
Also, the sharpness of the peak at the Fermi energy re- 
markably changes with the concentration. 

The main features of the DOS profiles for ordered com- 
pounds are assigned by Zou and Fu,lI in the order of in- 
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FIG. 8. Theoretical KKR-CPA (open circles) and exper- 
imental (dots) concentration dependence of the lattice con- 
stant and the magnetic moment in the Fe-Al alloys. 



creasing of energy, to the bonding, nonbonding (the peak 
just below Ep in NiAl) and antibonding states (Fig. 1^). 
A 'pseudogap' separating the bonding and nonbonding 
states remains well defined in all DOS calculations for 
ordered compounds. As is shown in Fig. S, the non- 
bonding structure completely disappears for disordered 
alloys. The tendency for the disappearance of the non- 
bonding peak can be already traced in Figs, g and || for 
the AS local DOS, but only the complete disorder leads 
to the total smearing out of this DOS feature. This effect 
is favorable for the magnetism, which develops in Co-Al 
disordered alloys. On the other hand, the disorder re- 
duces the formation energy (see Table O) because the Al 
to TM charge transfer requires that a late TM atom has 
as many Al atoms as its nearest neighbors as possible, 
in order to facilitate the charge transfer and the bonding 
hybridization. 
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FIG. 9. Total DOS for several concentrations of disordered 
Fe-Al, Co-Al and Ni-Al alloys as calculated by KKR-CPA. 

over the range of compositions 75 at.% to 40 at.% of Fe, 
i.e. the existence range of the B2 phase. In Co-Al dis- 
ordered alloys, the onset of ferromagnetism occurs only 
for a Co concentration larger than 50 at.%. The Ni-Al 
alloys do not show any magnetic ordering. 

The Friedel oscillations around the defects are favor- 
able for the enhancement of magnetic moments in the fer- 
romagnetic FeAl matrix. It means that any defect creates 
an enhancement of the magnetic moment on neighboring 
Fe atoms, as charge transfer to these atoms decreases. In 
the paramagnetic Co-Al matrix, however, these oscilla- 
tions are only sufficient for creating magnetic moments 
on Co sites, which are directly neighbors to the Co AS 
impurities. 



ACKNOWLEDGMENTS 

This work was supported by the Deutsche Forschungs- 
gemeinschaft. N.I.K. thanks the University of Osnabriick 
for the kind hospitality during his stay there. Useful com- 
ments by J. Kudrnovsky are appreciated. 



V. CONCLUSION 

We have shown in this paper that the onset of mag- 
netism in the Ni, Co and Fe aluminides is closely re- 
lated to the defect structure of these compounds. Among 
the perfectly ordered 1:1 compounds, only FeAl retains a 
magnetic moment of 0.76 fig per atom. In CoAl and NiAl 
the magnetic moments are totally quenched due to the 
d-band population. However, in the case of TM-AS de- 
fects the value of the magnetic moment for a Fe defect is 
the same as for a pure bcc iron atom. A large magnetic 
moment also appears for the Co- AS defect. All nearly 
disordered Fe-Al alloys exhibit a ferromagnetic ordering 
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